Numerical study of mesoscale cellular convection

Mesoscale cellular convections over the East China Sea during cold air outbreaks are simulated with a high-resolution numerical model. The model incorporates important physical processes involved in shallow convection, such as the exchange of heat and moisture between water and air; condensation; evaporation; and vertical turbulent transfer of heat, moisture, and momentum.The results show that open cells develop with aspect ratios as large as 14. The structure of the convection is examined in detail. The organized mesoscale circulation is responsible for breaking up the initial stratus cloud deck and enhancing turbulence in the upward-moving area (especially inside cumulus clouds). However, it is found that the heat flux contributed by MCC's themselves is much smaller than the eddy heat flux.     

Correlation of foreshocks and aftershocks and asperities

A close correlation in spatial distribution of local seismic activity and energy release patterns before and after the 1979 Petatlan, Mexico earthquake suggests heterogeneity within the fault plane of this major low-angle thrust event associated with subduction along the Middle America Trench. A simple two-asperity model is proposed to account for the complexity. Foreshocks and aftershocks of the neighboring 1981 Playa Azul earthquake showed a similar pattern. As both events occurred at the junction of the Orozco Fracture Zone and the Middle America Trench, we speculate that the observed complex fault plane is caused by subduction of the rugged ocean floor of the Orozco Fracture Zone. Short-term precursory seismicity prior to the Petatlan earthquake can be explained by using the asperity model and migration of a slip front from the south-east to the north-west across the main shock source region.     

Effect of cold-air advection on internal boundary-layer development over warm oceanic currents

Field experiments have shown that an internal boundary layer (IBL) can develop when the wind blows perpendicularly from the cold-water side of a warm oceanic current to the warmer side. The height of this IBL downwind from the thermal front as computed by an equation originally developed by Venkatram is in excellent agreement with observations.     

A derivative-based interpretation approach to estimating source parameters of simple 2D magnetic sources from Euler deconvolution, the analytic-signal method and analytical expressions of the anomalies

The major advantage of using either the analytic‐signal or the Euler‐deconvolution technique is that we can determine magnetic‐source locations and depths independently of the ambient earth magnetic parameters. In this study, we propose adopting a joint analysis of the analytic signal and Euler deconvolution to estimate the parameters of 2D magnetic sources. The results can avoid solution bias from an inappropriate magnetic datum level and can determine the horizontal locations, depths, structural types (indices), magnetization contrasts and/or structural dips. We have demonstrated the feasibility of the proposed method on 2D synthetic models, such as magnetic contacts (faults), thin dikes and cylinders. However, the method fails to solve the parameters of magnetic sources if there is severe interference between the anomalies of two adjacent magnetic sources.     

Prediction of surface horizontal displacements, and gravity and tilt changes caused by filling the Three Gorges Reservoir

 Horizontal displacements, and gravity and tilt changes induced by filling the Three Gorges Reservoir are modeled using elastic loading Green functions. When the water surface reaches its highest level, the effects become maximum on the reservoir banks. The longitudinal and latitudinal components of the horizontal displacements reach −8.2 and 7.7 mm respectively, gravity is increased by up to 3.4 mGal, and the prime vertical and meridian components of the tilt changes are −7.8 and −17.5 arcseconds respectively. Accordingly, the filling of the reservoir will influence values observed from global positioning system (GPS), gravimetry and tilt measurements in the area. The results given can be used to provide important corrections for extracting earthquake-related signals from observed data. Received: 19 January 2001 / Accepted: 3 September 2001     

The stability of merwinite in the system CaO-MgO-SiO2-H2O-CO2 with CO2-poor fluids

The stability of merwinite (Mw) and its equivalent assemblages, akermanite (Ak)+calcite (Cc), diopside (Di)+calcite, and wollastonite (Wo)+monticellite (Mc)+calcite was determined at T=500–900° C and P _f=0.5–2.0 kbar under H₂O–CO₂ fluid conditions with X _{CO 2}=0.5, 0.1, 0.05, and 0.02. Merwinite is stable at P _f=0.5 kbar with T>700° C and X _{CO 2}<0.1. At P _f=2.0 kbar, the assemblage Di+Cc replaces merwinite at all T and X _{CO 2} conditions. At intermediate P _f=1 kbar, the assemblage Ak+Cc is stable above 707° C and Wo+Mc+Cc is stable below 707° C. The univariant curve for the reaction Di+Cc=Wo+Mc+CO₂ is almost parallel to the T axis and shifts to low P _f with increasing X _{CO 2}, with the assemblage Di+Cc on the high-P _f side. The implications of the experimental results in regard to contact metamorphism of limestone are discussed using the aureole at Crestmore, California as an example.     

LSQRD及解的分辨分析

较详细地介绍了ＬＳＱＲＤ（Ｌｅａｓｔ Ｓｑｕａｒｅ ＱＲ Ｄｅｃｏｍｐｏｓｉｔｉｏｎ） 反演方法以及“检测版”分辨分析方法,并举实例检验其应用。对于大型稀疏方程组的求解,ＬＳＱＲＤ 是一种内存需求小、计算效率高和抗噪能力较强的算法。“检测版”试验表明,ＬＳＱＲＤ也可以对迭代解进行有效的可靠性估计,以指导对反演结果的利用。     

Wavelet filter analysis of atmospheric pressure effects in the long-period seismic mode band

The importance of the reduction of atmospheric pressure effects becomes very clear when investigating seismic normal-mode spectra below 1.5 mHz. The usual simple correction method consists in subtracting a term converted from local atmospheric pressure (pressure multiplied by a frequency-independent admittance) from the gravity record in time domain. Thus, estimating an efficient admittance is the key for an improved correction. Band-pass filters derived from dyadic orthogonal wavelet transform, having narrow pass-bands with good frequency response but without Gibbs phenomenon and causing no phase lag, are very helpful to estimate an efficient admittance, which is both time and frequency-dependent. Processing of high quality superconducting gravimeter (SG) records for the great Sumatra earthquake (Mw = 9.3, Dec 26, 2004) with wavelet filters reveal the three very well resolved splitting singlets of overtone ₂S₁ with a single gravity record after correction with time-dependent and frequency-dependent admittances. We also observe all coupled toroidal modes below 1.5 mHz, except ₀T₅, ₀T₇ and ₁T₁, with good signal-to-noise ratio (SNR); moreover, toroidal modes ₁T₂ and ₁T₃ are for the first time unambiguously revealed in vertical components.     

Application of statistical method in pollution episode analyses

Observed pollutants are of both primary and secondary origins, influenced by local topography and meteorology. Identifying sources and relative contribution often require vast detailed data collection and complicated models. This study applied a statistical time series analysis to two selected pollution events, spring and fall, at two sites in northern Taiwan. Vector moving average representations were used to quantitatively examine relationships among chemical pollutants and estimate their lifetimes. Results from impulse responses show that wind direction change alters the characteristic of pollution observed in opposite sites of Taipei City, from chemical dominant system to transport dominant one and vice versa. Chemicals are clearly separated into photochemical pollutants and primary pollutants. Results pointed out that Taipei City is the major source of photochemical smog, but not these primary pollutants. Derived chemical lifetimes at same location vary from 20% to four times under different meteorological condition. Estimated concentrations of hydroxyl radical range between 2 to 8 × 10⁷ cm⁻³. Photochemical pollutants are responsible for parts of PM10 collected in both station observed. Oxidation of SO₂ is only important in PM10 observed at one station. This study provides a simpler tool to derive information usually from complex models, therefore, is suitable as complement in decision-making process.     

水准椭球的纬向密度分布

本文推导了旋转椭球极点重力与赤道重力的密度积分公式。按照１９８０大地参考系统及水准椭球的“极点重力条件”和“赤道重力条件”,从数学上证明了地球的“赤道径向平均密度”大于“极点径向平均密度”,并给出了水准椭球纬向密度分布的估算公式;进而提出了“地球的径向平均密度按纬向正常分布”的假说。